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Abstract: In this paper, we shall address some field theoretic issues regarding the chiral magnetic 
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careful examination of the infrared limit of the relevant thermal diagrams. 
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1. Introduction 

The chiral magnetic effect (CME) proposed in ^, ^ provides a new probe of the QCD phase 
transition and the formation of quark-gluon plasma via relativistic heavy ion collisions(RHIC). The 
physical picture of CME relies on the interplay between the helicity of a quark and the external 
magnetic field. Consider a quark with positive (negative) helicity, its magnetic moment and the 
electric current it carries are always parallel (antiparallel) independently of the sign of its electric 
charge. The magnetic moment tends to be parallel to the magnetic field, so the electric current 
will be parallel (antiparallel) to the field for positive (negative) helicity. For massless quarks, the 
helicity coincides with the axial charge, 

Q5 = y d'r^7475V' (1-1) 

with the quark spinor tp carrying both color and flavor indexes. Therefore, for QGP of a nonzero 
axial charge density, a net electric current will be generated in (opposite to) the direction of the 
external magnetic field if the positive (negative) helicity is in excess. 

The conditions that support CME are likely implemented in RHIC. Firstly, for off-central 
collisions, a strong magnetic field is produced perpendicular to the collision plane; Secondly, because 
of the high temperature, there may be a sizable probability for the transition to a topologically 
nontrivial giuon configuration accompanied by a change of the axial charge according to the winding 
number ||, ||, |6| 

where Fj^^ is the strength of the color SU (Nc) field {Nc = 3) with / the color index and Nf is the 
number of flavors. Thirdly, the de-confined quarks that carry the chiral magnetic current can travel 
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sufficiently far before hadronization to lead to observable charge asymmetry perpendicular to the 
collision plane. It has been suggested recently that such a charge asymmetry is correlated with the 
baryon number asymmetry through a similar mechanism, the chiral vortical effect |0, ||] . For the 
experimental status of CME, see for example g |l^, |n|, |l| 

The chiral magnetic effect for a free quark gas in a static and homogeneous magnetic field B 
at thermal equilibrium has been analyzed in great details. With the aid of the grand partition 
function at a nonzero axial chemical potential /is, 

Z = Tre-'^(^-^^-''='5=) (1.3) 

with H the Hamiltonian, N the quark number, (3 the inverse temperature and fi the quark number 
chemical potential, one obtains the chiral magnetic current 3 — r]j where 

^ = ^cE'?/ (1-4) 
/ 

with Qf the charge number of the flavor / and the current per unit charge given by the classical 
expression 

j = |^M5i3. (1.5) 

The chiral magnetic current at nonzero momentum and frequency has also been calculated via 
current-current correlator to one loop order within the same grand canonical ensemble defined by 
( O ) Q . The same effect has also been examined with holographic models [|l^ |l^, |l^, |l^, |l^, |l^, 



and the lattice simulation pu, E2l . The effect of a nonzero quark mass has been considered recently 



in [g3|. A diagrammatic proof of (1.5) to all orders at high density has been attempted in ]24[ ]. 

It was pointed out in ||l^ that the naive axial charge ( p^.l| ) is not the right object to define the 
grand canonical ensemble since it is not conserved because of the axial anomaly. 



dx^ ' 327r2 ''^-P^^t^-^P^ + *^i67r2 Wa^I^.^^pA - , (1.6) 

where the axial vector current Js^ = i^pj^j^ip ^-^d is a linear combination of the Chern-Simons 
of QCD and QED, given by 

n - .A' f^-l f^'^'^A'^A + ^v—e .A ^ fl 7) 

8^2 WA^. (^^^^ 3/ A^A^j+ifj^^^e^^.^A^g^^. (1.7) 



with A and A^ the gauge potential of gluons and photons. The integration of (1.6) gives rise to 



(1.2), to which the trivial topology of the electromagnetic field does not contribute. The conserved 



axial charge to replace Qs in (1.3) reads 

(35 = (35+* / rf^rf^4. (1.8) 



In what follows, we shall name Qr, the naive axial charge. Furthermore, the author of argued 
that the gauge invariance prevents a nonzero chiral magnetic current to be generated from the grand 
canonical ensemble defined with and the chiral magnetic current comes solely from the second 



term of (1.8) in the ensemble defined by Qq. Because this term stems from the anomaly, which is 



universal to all orders, the classical expression (1.5) is robust against higher order corrections 



In this paper, we shall analyze the chiral magnetic effect via the current-current correlator 
in the light of Rcf. ||l^. There are standard recipes to implement gauge invariant regularization 
schemes via thermal diagrams employed in this work. Higher order corrections can also be included 
systematically. We find that the validity of the statement in relies on the existence of the 
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infrared limits of the energies and momenta involved, which is not always guaranteed. We shall 
pinpoint a few exceptions to the statement in jl^, one is caused by the massless poles of the 
invariant form factors underlying the triangle diagram at T = and /i = and others are related to 
the noncommutativity between the zero momentum limit and the zero energy limit at T 7^ and/or 
/Lt ^ 0. The latter subtlety is a common feature of thermal field theories. The difference between 
different orders of limits is likely to be subject to higher order corrections. Since the magnetic field 
in RHIC is neither homogeneous nor static and the system is not in a complete thermal equilibrium, 
these issues need to be addressed to assess the robustness of the effect in RHIC phenomenology. 

In the next section, we shall work out the most general structure of the chiral magnetic current 
consistent with the rotation symmetry, Bose symmetry and the gauge invariance. We shall restrict 
our attention to the diagrams that contribute to the same powers of /is and B as (L5). The 
infrared subtlety is allocated to some invariant form factors of three point functions. The one-loop 
evaluation of the chiral magnetic current will be revisited in the section III with the Pauli-Villars 
regularization. In the section IV, we shall clarify the relation between the chiral magnetic current 
and the axial anomaly for an inhomogeneous and time-dependent /is, which is related to the QGP 
off thermal equilibrium. The section V will conclude the paper with some open questions. 

Throughout the paper, all four momenta will be denoted by capital letters. We shall adopt the 
Euclidean metric (1, 1, 1, 1) in which a Minkowski four momentum P = (p, ipo) with po real. All 
gamma matrices are hermitian. 



2. The General Structure of the Chiral Magnetic Current 

The Lagrange density of a quark matter at nonzero baryon number and axial charge densities is 
given by 

^ = -IKA- - 7^M-^M- -^(^^lir- ^9T'A\ - zeqA^ ^ (2.1) 



+ /iV^74?A + /is (V'7475V' + i^i) + 

+ gauge fixing terms and renormalization counter terms 

where q is the diagonal matrix of electric charge in flavor space, /i is the quark number chemical 
potential and /is is the axial charge chemical potential. An external electric current J^^*"- has been 
added to the Lagrange. 

The generating functional of the connected Green function of photons is the logarithm of the 
partition function 

Zir""^-] ^ [ [dA^][dA][d^][dip]exp (i [ dtd^rc] . (2.2) 



For the Matsubara Green functions, the time integral inside exp(...) is along the imaginary axis 
of the complex i-plane extending from to i/3 ~ i/T subject to periodic (antiperiodic) boundary 
conditions for bosonic (fermionic) fleld variables and —T\nZ is the thermodynamic potential at 
equilibrium. For the closed time path Green function (CTP), the time t is integrated along the real 
axis from — cxd to 00 and then from 00 back to — c» and the thermal equilibrium is implemented by 
the initial correlations. All flelds can take values on either branch of this contour, which doubles the 
number of degrees of freedom |2^, ^ ^ . See appendix ^ for a brief introduction of the CTP 

formalism. The external current J^^^' generates a nonzero thermal average of the electromagnetic 
potential, given by 

. / N SlnZ ,^ 
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and its Legrendre transformation reads 



where the effective action 
S[A] 



ilnZ[J'=^'- 



(2.4) 



(2.5) 



7^ ilV-F {ly 



with Tft,v = — and Bi — ^eijkJ'jk = (V x A)i. In the second hne of (2^), we have 
separated the contributions of tree diagrams (first two terms) from that of loop diagrams (third 
term). Eq.(2.4) is equivalent to the Maxwell equation 



J?. 



where 



Ji{x) = 



ST 



Jut 



(2.6) 



(2.7) 



6Ar{x) 

represents the induced current in the medium. The functional T[A] can be expanded according to 
the powers of A with the proper vertex functions as coefficients. We have, in momentum space 

d^Q r 1 



T[A] 



(27r)4 



-n^,(Q)^;(Q)A(g) + o(A) 



(2. 



where only the term contributing to the linear response is displayed explicitly. It follows from (2.7) 
that 

J.(0) =/C„(Q)A(Q), (2.9) 



where 



ICij{Q) = -Ily (Q) - irj^-^^isCijkqk + 0{A^ 



(2.10) 



with all QCD and higher order QED corrections contained in the photon self-energy tensor I\-^u{Q)- 
The prescription of the functional derivative for the retarded linear response is outlined near the 
end of the appendix The antisymmetric part of ICij{Q), 



(2.11) 



which is odd in n^, carries odd parity and generates the chiral magnetic current. 

'-fj iQ) the powers of /Z5 , we have /C^ 

=2 



Expanding the response function ICfAQ) in the powers of we have ICfAQ) — A*5^.|}''(<3) 



0(jU5), where 



(2.12) 



and underlies the classical form of the chiral magnetic current (1.5). 

The first term of (2.12) is represented by the IPI diagram with two external vector vertices 
and an external axial vector vertex, shown in Fig.l, at 11 — i, — j and p = 4. The lowest order 
of it consists of the usual triangle diagrams in Fig. 2. Let the incoming 4-momcnta at the photon 
vertices be Qi = (qi,ia;) and Q2 = (q2,— iw), the incoming 4-momentum at the axial vertex is 
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Figure 1: The diagrammatic representation of the contribution to the chiral magnetic current from the 
photon self-energy, where the contribution of each vertex to the Feynman amplitude is indicated explicitly. 




then — Qi — Q2 = (^qi — q2, 0). The amphtude of the diagram A^i,(Qi, (52) consists of a pseudo- 
tensor Aij((3i, Q2), a pseudo-vector A4j{Qi,Q2) and a pseudo-scalar A44((5i,Q2)- In the limit 
Qi — >■ —Q2 with Qi = Q = (q, icj), we find that 

^n^.(g)|^,=o = Am.(Q, -Q). (2T3) 

The rotation invariance and the Bose symmetry 

A^.(Qi,g2) = A,^(Q2,Qi) (2.14) 
dictates the following most general tensorial structure 

Aij((3i,Q2) = ir?— [Co(gi,(7|,qi • q2;w)ejjfe(7ife - Co(g|,gi,qi • q2; -w)ejjfe(j2fc (2.15) 

+ Ci{ql,ql,(ii ■ q2; w)ejfci9ifeg2;9ii - C'i(g2,<??,qi ■ q2; -^^)<^ikiqikq2iq2j], 
A4fc(Qi,Q2) = il-:r^C2{ql,q2,<ii ■ (i2\^)<^ijkqiiq2] = ^k4{Q2,Qi) (2-16) 

Ztt 

and A44((5i, Q2) = 0, where Cq, Ci and C2 are dynamical form factors. The time reversal invariance 
implies that Cq, Ci are even functions of w and C2 is odd in w (This, however, is not required for 
our purpose). Notice that the tensors ^ikiqikq2iq2i and eikiqikq2iqij are not independent and can 
be reduced to the tensors already included in ( 2.15 ) via Schouten identity 

^ijkqi — ^lijqk + ^knqj — ^jkiqi = 0. (2-17) 
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Furthermore, switching Qi — > Q2 amounts to qi — S- q2 and lo —oj. It follows from (2.12) and 



(12.151) that 



l^fjiQ) = ^V^f^5[F{Q) ~ l]e,jkqk + O(pI) (2.18) 

with 

F{Q) = -Co(g^ ~q';Lu) - C^{q\q\ -q^; -u:). (2.19) 

The chiral magnetic current in a constant magnetic field corresponds to the limit -F'(O), which is 
subtle as we shall see. 

The electromagnetic gauge invariance, 

QipA^,(Qi, Q2) = g2.A^,(Qi, ga) = (2.20) 

gives rise to the relations 

Co{q1,q2,<li ■ q2;w) = -qlCi{ql,q1,Cii ■ q2; -w) + a;C2(gi, q?, qi • Qa; (2.21) 

and 

Ca{ql,ql,qi-q2;-i^) = -qlCi{ql, qj, qi ■ q2; uj) - ujC2{ql, qj, m ■ q2; (2-22) 
and therefore 

F{Q) = q^[C^{q\ q\ ~q^; Oj) + C,{q\ q\ -q^; -uj)] + c.[C2(q^ q\ -q^; c.) - C2{q\ q\ -<z'; 

(2.23) 

If the infrared limit of the dynamical form factors C\ and C2 exists, then i^(0)=0 and there 
is no chiral magnetic current associated to the naive axial charge. This is the case in the static 
limit q — > with Q = (q, 0) to one- loop order at nonzero T and/or /x. It remains so if there exists 



an nonperturbative IR cutoff to remove the \ singularities brought about by QCD corrections! 30 



(Such kind of singularities is likely to occur for diagrams with more than one quark loops linked by 



gluon lines). In that case, the chiral magnetic current takes the classical form (1^) to all orders. 

It is a common feature of thermal field theories that the different orders of the double limits 
limq_s.o lini^_>o and limtj_>.o limq_j.o may not agree. While the former order of limits of Ci(g^, g^, — (j^,w) 
and C2(g^,(7^, — g^; w) converges and leads to the classical form of the chiral magnetic current, the 
latter order of limits leads to IR divergence. The explicit calculation of the triangle diagram of 
Fig. 2 in the appendix ^ with /i = 4, p = 4 and v ^ j yields 

C2(0,0,0;l.) = -^ (2.24) 

as — J> and limtj_j.o lim^-s-o ^(Q) = §• Consequently, the magnitude of the one- loop chiral 
magnetic current is reduced to one third of the classical magnitude. This is consistent with the 
direct one-loop calculation in the literature [T^ and will be reexamined in the next section. Since 
the form factor F{Q) is not linked to the axial anomaly, the chiral magnetic current in this order 
of limits is likely to be subject to higher order corrections. 

The IR singularity also shows up via the massless poles if the zero temperature and zero chemical 
potential limits are taken prior to the limit Q ^ and A^,i,(Qi, Q2) becomes fully covariant then. 
To the one-loop order, the triangle diagram Fig. 2 gives rise to 

gi(g'.g^V;^)^ ^(^2^^2) (2-25) 

and 

C2{q\q^-q';L0)^~ (2.26) 

(See section IV for details.) Both Ci and C2 are infrared divergent and we find F{0) = 1 and 
therefore zero chiral magnetic current for T = /i = but /is 7^ 0. 
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p 

Figure 3: The one-loop diagram of the photon self-energy. The solid line with a double arrow stands for 
the free propagator to all orders of 



3. The one-loop contribution 

The one-loop contribution to the chiral magnetic current has been discussed extensively in the liter- 
ature. In the present section, we shall supplement this calculation with the Pauli-Villars regulariza- 
tion, since the photon self-energy as a whole suffers from the UV divergence. As the regularization 
respects the gauge invariance, the result will be consistent with the Ref.[|l^ and the statement of 
the previous section. The trivial color-flavor factor rj will be suppressed below. 

The one-loop photon self-energy tensor at the temperature T, shown in Fig. 3 is given by 

U^AQ)^e^Tj2 [ T^[s^.(P,Q|m)-^C,S^,(P,Q|M,)], (3.1) 

where 

S^.(P,Q|™) = tTSF{P + Q\m)j^SF{P\m)j,. (3.2) 

and the summation in the integrand corresponds to the contribution of the Pauli-Villars regulators 
that remove all UV divergences. We have 

Y^Cs^l (3.3) 

s 

and Ms — > oo after the integration. The free quark propagator with a four momentum P = (p, ipo), 
a mass m, a quark number chemical potential fj. and an axial charge chemical potential /i5 reads 

SAP\m) = — (3.4) 

/P + Ai74 + Ai57475 - m 

= -^[MPim, fi, H5) + MP^m, H, -fJ.5)] + ■^l5[B{P,m,n,fi5) - B{P,m, n, -^5)] 
where /P = 74P0 — 17 • p and we have decomposed S'i?(P|m) into the parts even and odd in /15 with 

A{P, m, fi, /is) = — -—2 ^— ^2 — (3-5) 



and 

-{ p + ^5)74 + »(po + ^J■- ^74) 7 • P 

{po + /i)2 - {p + ^izY - m 
The chiral magnetic current corresponds to the antisymmetric spatial components of Ilf^^{Q), i.e. 



B(P, m, 11, H5) = , ,,,2 _ / , X2Z — 2 ■ (3-6) 



n^(Q)EEl[n,,(Q)-n,,(g)] = e2r^ / ^[s,^.(F,Q|m) -5]c,s^(f,q|m,)], (3.7) 
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where 

S^(P,Q|m) = ^^tv-f5{[BiP + Q,m,^i,fi5)-B{P + Q,m,fi,~fi5)h, (3.8) 

+ [B{P, TO, Ms) - B{P, TO, fi, -fi5)]'-fi [A{P + Q, TO, fi, Ms) + A{P + Q, to, fi, -fj,5)hj}- 

It is straightforward to work out the trace and summation over the Matsubara frequency, po = 
i{2n + l)7rT. To obtain the retarded self-energy, we shall follow the recipe of Baym and Mermin 
||3l||to extend the Matsubara frequency qg to the upper edge of the real axis, qo ^ uj + iO+. The 
details are shown in the appendix ^ and we shall report two special cases below. The antisymmetric 
part of the self-energy tensor is parametrized as 

= -«^A^5Fi(g,^)e*,fcgfc, (3.9) 

with Fi{q,Lu) at fi^ — corresponds to the one-loop approximation of F{q,uj) as defined in Eq. 
( |2.19|) . The dependences on the spatial momentum and the energy are indicated separately here. 
Diagrammatically, Fig. 2 corresponds to the linear term of the Taylor expansion of Fig. 3 in fi^. 

3.1 The static limit 

At zero frequency, <?o = 0, we find that 

Fi(g,0) = -nq\m) +J2CsHq\Ms) (3.10) 

s 

where 



J'(glTO) = 



1 



dppln 



1 

'2p + gl 



(3.11) 



{P±J^[f{E+ -^^)- f{-E+ m)] - ^^[fiE^ M) - f{-E^ m)]} 



with 



E± = ^{p ± Ais)' + m2 



and the Fermi distribution function 



/(0 = 



1 



6/9? + 1 ■ 

It is straightforward to verify that the limit g — > at T 7^ and/or /i 7^ yields. 



(3.12) 
(3.13) 



HO\m) I dp{[ 



2/^5 Jo 

1 



, P + 

E+ 



[f{E+ -^i)- f{-E+ - ^^)] - 



E- 



[/(i?_-M)-/(-i?--M)]} 



2/i5/3 
1 



ln(l + e-'3(iJ+-f)) _ in(i + e-/3(S— m)) + ln(e'3(^++") + 1) - ln(e^(^-+^) + 1) 



lim(ii;+-i?_) = l, 

Z/i5 p-i-oo 

Then eqs.( |3.1C| ) and (3^) implies that 



lim lini Fi (g, uj) = 0. 



(3.14) 



(3.15) 



This result is expected according to the discussion in the last section because the nonzero Matsubara 
frequency, {2n + l)TrT regularizes the infrared behavior of the quark propagator even in the massless 
limit. Notice that the regulator contribution 



M. 



lim ^(g|Af,) = l 



(3.16) 
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for all q and this is also the case with a non static Q. If, on the other hand, T and /i as well as the 
quark mass are set to zero first, we find 



Ai5<? Jo 



dppln 



2p-q 



(3.17) 



It follows from ( |3.16D , ( |3.17D and ( |3.10D that Fi{q,0) 
covariant result reported at the end of the last section. 

3.2 Massless limit 



2p + q 

— 1 at = 0, in agreement with the 



In the massless limit, m = 0, the quark propagator (3.5) reduces to 
Sf{P\0) = 



1 + 75 



1-75 



and Fi(q,uj) in this case reads 



Fi{q,Lj) 



2fi. 



5 Jo 



dpp' 



J {p, q,^^) + J {p, -q, -w) J (p, 1, -w) + J{p,-q,uj) 



(3.18) 



e/3(p+p-P5) _|_ 1 



J{p,q,^) + J{p,-q^-^) , J(p,q,-^) + J{p,-q:^) 



1, 



where 



and 



ReJ{p,q,uj) 



pq 



UJ 1 

1^2 



1+UJ 



Im J(p, g, uj) ~ — sign(w) 1 + cj 



— 2puj - 
2p(f- 

LtP' — 2puj — q^ 



In 



(g; - q){uj + q - 2p) ^ 
{uj + q){oj - q - 2p) 

\q^ + 2puj — u}^\ 



1 - 



2pq 



(3.19) 



(3.20) 



(3.21) 



pq \ 2pq^ 

and the "+1" of cq.( ^.19 ) comes from the Pauli-Villars regulators. The limit Q — s> of the PV 
regulators is independent of the order between q — > and a; — > as long as Ms — > cxd is taken first. 
The same limit of massless part is, however, subtle. We have 



but 



lim lim Fi(q,uj) — 

g->Oaj->0 



lim limFi(g,a;) — — 

cx;->-Og->0 3 



(3.22) 



(3.23) 



consistent with the result reported in The nonzero value of the latter limit signals infrared 

divergence of the form factor 6*2(9^, g^, —q^\ w) defined in the last section under the same orders of 
limits. 



4. The Relation to the Triangle Anomaly 

In the section 2, we related the chiral magnetic current to the infrared limit of the three point 
Green's function in Fig.l with two electric currents and the fourth component of the axial vector 
current. We analyzed the general structure of the chiral magnetic current as is required by the 
electromagnetic Ward identity. For the sake of simplicity, we restricted our attention to zero energy 
flow at the axial vector vertex. To explore the the impact of the anomalous axial current Ward 
identity, this restriction will be relaxed in the present section. The physics of the diagram of Fig. 1 
with p = 4 and an arbitrary Qi + Q2 corresponds to CME at a space-time dependent /i5 in a QGP 
off thermal equilibrium. 
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We shall denote the general Feynman amplitude of Fig.l by h^^upiQiTQi) with Qi and Q2 the 
incoming momenta at the vector vertices indexed by /i and v. We have 

A^.4(Qi,Q2) = A^,(Qi,Q2) (4.1) 

with Api/((3i,Q2) defined in the section 2. The incoming momentum at the axial vector vertex is 
then 

i^ = (k,zfco) = -Qi-Q2- (4.2) 

We have 

QimA^i.p(Qi, Q2) = Q2uKyp{Qi^Q^) = (4.3) 
following from the electromagnetic gauge invariance. The triangle anomaly implies that 

(Ql + Q2)pApyp(Qi, Q2) = -'ir]Tr^'^^iiyat3QlaQ2l3 (4.4) 



which holds to all orders of interaction at arbitrary temperature and chemical potential ||3^ . The 
classical expression of the chiral magnetic current is associated to the component Aij^lQi, Q2) with 
the momenta 

Qi = (q, iw) Q2 = (-q, -iuj). (4.5) 

It is tempting to relate the self-energy contribution to CME with the axial anomaly via the 
limiting process 

1 e2 
A.y4(Qi,Q2) = -i lim ^(g'l + <32)pAw(<9i,<32) = -iVTT^^ijkqk (4.6) 

where Q'l = (q, ifco/2) and Q'2 = (— q, «/(;o/2). This appears in contradiction with the statement 
of the absence of CME with the naive axial charge. It does not display the nontrivial energy- 
momentum dependence of the one-loop result. The reason lies in the infrared singularity and the 
subtlety of the order of limits feg — and k ^ as we shall analyze below. At T = and /i = 0, 



however, the order of limits is irrelevant and we always get the RHS of (4.6), consistent with the 
one loop result near the end of the subsection 3.1. 

The most general tensorial decomposition aX T — ^ ~ Q consistent with the gauge invariance 



(4.3) and Bose symmetry reads 



A^i.p(Qi, Q2) = i'n^^{e^^uai3QiaQ2f}[QipDi{Q\,Ql,Qi ■ Q2) + Q2pDi{Ql,Ql, Qi ■ Q2)] 

+ {^vpal3QlaQ2pQlfj. — Ql£tJ.up\Q2x)D2{Qi,Q2TQl ■ Q2) (4.7) 
— {^Hpal3QlaQ2l3Q2y — Q2^piypxQl\)D2{Ql, Qi, Ql ■ Q2)}, 

where the 4D Schouten identity 

^^vpxQoL ^afiupQ X ~t- f^Xa^uQ p ~t- ^pXapQv ^vpXaQp — (4.8) 

is employed to reduce the number of terms. It follows from the anomaly equation ( |4.4| ) that 

(Ql + Q2) ■ QMQlQlQi ■ Q2) + (Ql + Q2) • Q2A(Q2, Ql, Ql • Q2) 
- QlD2{Ql,Ql Ql • Q2) - Q^^2(Q^, Ql Qi ■ Q2) = -i, (4.9) 

which implies infrared singularities of the dynamical form factors Di and D2. To the one loop 
order, we find that 
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Figure 4: The CTP diagram with one vertex insertion highlighted. 



and 



dx 







dy 



x{l - x-y) 







Qjx + Qly-iQ.x-Q^y)^' 



(4.11) 



which satisfy the constraint (|4.9|). For the CME momenta, ([l.5|), we find that D^iCt ,Q\ ~Q ) 
2^ and therefore 



e 



(4.12) 



Breaking the tensor (4.8) into spatial and temporal components, we obtain (2.25) and (2.26|) via 



<M3 



At a nonzero temperature and/or chemical potential, the limit K ^ becomes very subtle. 
Because of the discreteness of the energy in the Matsubara Green's function, one has to switch to 
the real time formalism for the analysis, of which, the closed time path (CTP) Green's function 
is most convenient. The main ingredients of CTP is summarized in the appendix ^ Explicit 
calculations of the triangle diagram via the CTP show that 



lim lim Aij4{Qi,Q2) ^ lim lim Aj.4((5i, (52). 



(4.13) 



with k and kg defined in (4^). The limit order on RHS leads to (4^), the result dictated by 
the anomaly, while the limit order on LHS gives rise to result of the last section, obtained from 
the Matsubara formulation and its analytic continuation to real energy. Therefore, there is no 
contradiction between the universality of the anomaly and the statement of ||l^ . 

The subtlety of this infrared limit can be explored in general. Consider the CTP diagram 
in Fig. 4 with a vertex insertion of four momentum K — (k, z/cq), summing up both CTP paths. 
The amputated external legs pertaining to the shaded bubble are suppressed. It follows from the 
Feynnian rules of CTP that the contribution of the two highlighted lines adjacent to the vertex 
insertion in the Fig. 4 is 



SiaiP + K)TSbl{P) - S2a{P + K)TSt,2{P), 



(4.14) 



where Sah{P) is the CTP quark propagator defined in the appendix ^ with a, b labeling the two 
CTP paths and F is a matrix with respect to the spinor indexes. The spinor indexes as well as the 



indexes a and b of ( 4.14 ) are to be contracted with the contribution from the shaded bubble of Fig. 
4. In terms of the retarded(advanced) propagator S r{P){S a{P)) and the correlator Sc{P) defined 
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in Eq.( 



we find that 



Sia{P + K)rSbl{P) ~ S2a{P + K)TSb2{P) 



(4.15) 



= -[Sc{P + K)rSRiP) + SAiP + K)rSc{P) ± SniP + K)TSr{P) ± SAiP + K)rSA{P)]. 

with "±" on RHS depending on the CTP indexes, a and b. Therefore the amplitude of the diagram 
has the foUowing mathematical structure 



g{K) 



(2^)4 

y(po,p;fco,k)} 



{c^(po,p;fco,k) 



[l~2f{p'm(P + Kf] , [1 ~ 2/(p)]<5(p2) 



p2 



{P + Kf 



(4.16) 



where = |p + k| and f{p) stands for the fermion distribution function. For the sake of simplicity, 
we have set the quark number chemical potential /z = 0, but the generalization to a nonzero ^ is 
straightforward. The quantity inside the bracket on RHS of (^ comes from the first two terms in 
the second line of eq.( 4.15| ) and the contribution from the shaded bubble of Fig. 4 is included in 
the functions U{pq, p; fco, k) and V{pq, p; fco: k). The function U {pa, p; fcoj k) is regular at the mass 
shells 

^'-p'-Po = {P + Kf =p'^ -p'o^ =0, (4.17) 

and its derivative with respect to pa will be denoted by C/(po, p; fco, k) below. So is the function 
V{pq, p; fcoi k) and its integral, I{K) is unambiguous in the limit K ^ 0. Carrying out the energy 
integral, we find that 



d^p 
(2^ 



{ 



1 



p-p' + ko 
1 



p-p' -k^ 



[1 - 2f{p')]U{p' - fco, p; fco, k) [1 - 2f{p)]U{p, p; fcp, k) 

p'{p + p' - fco) 
[l-2/(p')][/(-y-fco,P;fco,k) 



p{p + p' + fco) 



p'{p + p' + fco) 



[l-2/(p)]C/(-p,p;fco,k) 



} + I{K). 



It follows that 



lim lim QiK) 

/co->0k^O 



p{p + p' - fco) 

d^p l-2/(p)^C/(p-fco,p;fco,0) (7(p,p;fco,0) 



(4.18) 



- lim — 

2 feo-s-O fco 



(27r)3 



(2p - fco) 



(2p + fco) 



C^(-P- fco,P;fco,0) , C/(-p,p;fco,0) 



(2p + fco) 
d^p l-2/(p) 
(27r)3 p2 



(2p - fco) 
-J7(p,p;0,0) + 



f^(0)_ 

t/(p,p;0,0) 



C/(-p,p;0,0) 



/(O) 



and 



lim lim G{K) = — lim 



C/(-p,p;0,0) 

1 fl-2/(y)]t/(p',p;0,k) [l-2/(p)][/(p,p;0,k) 



(4.19) 



k^O fen-s-0 



2 k^o 7 (27r)3 p2 _p/ 



r{- 



[1 - 2fip')]Ui~p', p; 0, k) [1 - 2f{p)]U{~p, p; 0, k) 



p' 



P 



} + IiO) 



= lim lmig(K) + - ^ "^'^ d/C/(p,p;0,0) + C/(-p,p;0,0) 

fco^Ok^O 2 



(4.20) 



(27r)3 dp p2 

The inequality ( 4.13 ) is an example eq.( |4.2(]| ) for the three point function A.^i,p{Qi,Q2) with 



"*7574- Applying ( 1.20 ) for the one loop diagrams of Ay4(Qi,Q2) with ( |4.6| ) for the first 
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term in the third line, we recover the CME term of the photon self energy obtained previously. 
Indeed, the first term in the third line of ( 4.20| ) corresponds to the term on RHS of ( 3.19| ) 



and the integral of ( 4.20 ) in the same line goes to the integral of ( 3.19 ) in the limit /X5 — 0. Since 



the only /is dependence of (3.19) is through the distribution functions, the limit has its integrand 



proportional to the derivative of the distribution function. 
5. Discussions 

In this work, we investigated the interplay between the gauge invariance and the infrared limit 
in the chiral magnetic effect. The part of the induced electric current that is linear in the axial 
chemical potential /15 and the magnetic field B is divided into two terms, i.e. 

2 2 
HQ) = -Ti^ii5F{Q)B{Q) + ^^-^^,^B{Q) (5.1) 

where the first term corresponds to the loop diagrams of the photon self-energy tensor and the 
second term comes from the Chern-Simons term of the conserved axial charge Qs , which is dictated 
by the anomaly. The gauge invariance relates the form factor F{Q) to two form factors, Ci and 
C2 underlying a three point diagram of two vector current vertices and an axial current vertex. If 
the infrared limit of these form factors exists, F{Q) = to all orders of coupling and the classical 



form of the chiral magnetic current in a constant magnetic field, eq. (1.5) emerges. Our statements 
are illustrated with explicit one-loop calculations subject to the Pauli-Villars regularization. At 
zero temperature, however, both Ci and C2 are infrared divergent and i^(0) — 1. Consequently, 



the two terms on RHS of (5.1) cancel each other and the chiral magnetic current vanishes. At a 
nonzero temperature and/or a nonzero chemical potential, i^(0) depends on how the limit Q is 
approached. The magnitude of the chiral magnetic current is reduced if the zero momentum limit 
is taken prior to the zero energy limit, as is implied by the infrared divergence of C2 under the same 
order of limits. More subtle is the situation with a coordinate dependent /is. If the four momentum 
associated with /15, K = (k, ifco) is set to zero in the order limk-1.0 linifco-i-O: the results of sections 
2 and 3 are recovered. With the opposite order of the limit, however, F{Q) = 1 as is dictated by 
the anomaly and the two terms of ( ^.1| ) cancel again. Unlike what happens with the axial anomaly, 
the difference between different orders of the infrared limits is unlikely robust against higher order 
corrections. Since the ambiguity stems from quasi particle poles, it will disappear when the quasi 
particle weight is diminished by strong coupling. Then the chiral magnetic current will revert to 
its classical expression with the order livni^^o limg_i.o of the infrared limit Q ^ Q. This is consistent 
with the holographic result reported in fl^ . 

One complication with a coordinate dependent /is is that the term /isQs of the Lagrangian 



(2.2) is no longer gauge invariant. One may argue that this term is only defined in a specific gauge. 



say Coulomb gauge, in which the vector potential 



A = -i^VxB (5.2) 



is already gauge invariant. So is the Chern-Simons term of Q^. A possible objection to this approach 
is the violation of the micro causality, i.e. the commutator between two axial charge densities in 
Heisenberg representation does not vanish for a space-like separation because of the nonlocality 
introduced by the inverse Laplacian in (^|^). It remains an open issue to assess the validity of the 
conserved axial charge in a non equilibrium setup (See pq| for some related discussions). 



Finally, we would like to comment briefiy on the derivation of the classical result (1^) by 
summing up the single particle Landau orbitals in a constant magnetic field. It is a one-loop 
procedure to all orders of the magnetic field. The linear term of the electric and the magnetic 



- 13 - 



field stems from tlie same plioton self energy tensor discussed here and requires a gauge invariant 
regularization to cancel the UV divergence. In view of the analysis in this paper, we would expect 
that the summation over the Landau orbitals yields a null result for the chiral magnetic current if 
the regulator contribution is included. The net current is solely given by the Chern-Simons term of 
the ( |2.2D . Therefore we do not see any nontrivial effect of a nonzero quark mass claimed in . 
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A. Some elements of the closed time path Green functions 

The CTP formalism of a finite temperature field theory was introduced by Keldysh and 
Schwinger |2^. Good reviews can be found in The CTP contour on the complex 

time plane has two branches: Ci runs from negative infinity to positive infinity just above the real 
axis, and C2 runs back from positive infinity to negative infinity just below the real axis. All fields 
can take values on either branch of this contour, which results in a doubling in the number of 
degrees of freedom. The scalar propagator is given by, 

D{X -Y) = {T,<l>{XmY)) (A.l) 

where Tc is the operator that time orders along the CTP contour. We also use the notation 
X — (x, it) and P = (p, ipo). The propagator has 2^ = 4 components and can be written as a 2 x 2 
matrix 



D 

with 



Dii D12 
D21 D22 



Dn{x-Y) = (r(0(x)0(r))), 

Di2{X-Y) = {4>{Y)(^{X)) , 

D22{X -Y) = {fWX)(b{Y))), (A.2) 

where T is the usual time ordering operator, and T is the anti-chronological time ordering operator. 
These four components satisfy, 

Dii - D12 - D21 + D22 ^ 

as a consequence of the identity 6{x) + 9{—x) = 1. 

It is more useful to write the propagator in terms of the three functions 

Dr = Dn - D12 , 
Da = Dn - D21 , 

Dc = Dn + D22 ■ (A.3) 
Dn and Da are the usual retarded and advanced propagators, satisfying 

DrIX -Y)- Da{X -Y) = mX),<i>{Y)]) , 
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and Dc is the symmetric combination, also called correlator 

DciX-Y) = {{(j,iX),^{Y)}), 

which satisfies the Kubo-Martin-Schwinger (KMS) condition at thermal equilibrium. In momentum 
space 

i 

Dr,a[P) = 7 — I . — 2' 

Dc{P) ={l + 2n{po)){DR{P)-DA{P)) = 2n[l + 2n{Ep)]5{p-' + m'), (A.4) 
where n(po) is the thermal Bose-Einstein distribution, 

'^'^P^^ = efipo - 1 ' "(-Po) = - (1 + ri{po)) (A.5) 



and Ep — \/ p'^ + m^- The propagator can be rewritten asp£ 



2D = Dr (1, -1) + Da (_\) (1, 1) + Dc (I) (1, 1) . (A.6) 



Using the KMS condition (A.4) this expression can be rewritten as, 

D{p) ^ DR{p)(^\{l + n{p^),n{Pr>))-DA{p)(, J (1,1)- (A.7) 



1 + n(po) 

The fermion propagator S{P), can be obtained by multiplying D{P) with /P + m and replacing 
uIpq) with —fipo), the Fermi-Dirac distribution function 

We have 



S 



Su Si2 

S21 S22 



with Sii — S12 — 5*21 + S22 = 0. The retarded, advanced and correlator components are 

74Po -ij -p + m 



Sr = Sii — S12 — i 

Sa = Sii — S21 = i 



{po + «e)^ — — 
74P0 — i-f ■ p + ni 



{po — ie)^ — — ' 
5c = 5ii + S22 = 27r(74Po - ^7 ■ P + - 2f{Ep))5{P^ + m^) . (A.9) 

They satisfy the KMS condition 

Sc{P) ={l-2f{po)){SRiP)-SA{P)) (A.IO) 

We can extract the IPI two-point function, or self energy, by removing external legs in the usual 
way. We find, 

Ur = Hii + U12 
Ua = Hii + U21 

nc = nn + n22 (A.ii) 
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where II/j and Ha are the usual retarded and advanced self energies. The four CTP components 
satisfy the constraint, 



Hii + ni2 + n2i + n2 







The eq.(2^) in CTP formalism takes the form 



r[^] = 



\t\^.{qT''K*{q)aI{q) + o{a') 



5T 



U'=o= 



-fn" +ni2 _tt2i -n22w - _r[R a 



The recipe to obtain nonlinear responses is given in |2£ 



(A.12) 



where the superscripts a, 6, label the two CTP components. Introducing A — ^{Ai — A2) and 
A' ^ Ai + A2, we find that 



(A.13) 



B. The infrared behavior of the form factor C2('Zi,?2iQi ' ^2',^ 



In this appendix, we shall derive the infrared behavior ( 2.24[ ) by calculating the triangle diagrams 
in Fig.l with fi = p — 4 and v = j . The trivial color-flavor factor will be suppressed. We shall start 
with a nonzero Matsubara frequency lo = 2imTT and continue it to the upper edge of the real axis 
following Baym-Mermin prescription. The amplitude of the diagram reads 



A: 



1 



1 



1 



1 



1 



(B.l) 

where P = (p, i{po + /i)) and P± = (p. i(j>o + fi±uj)) and = —17 • q. Since we are only interested 



the form factor C2 at qi = — q2 = 0, we expand (BA) according to the powers of spatial momenta 
qi and q2 and pick up the term proportional to the product of them, i.e. 

A.^^N 2 f d^P ^ /^l/l 1 1/1 1/1 1 1/1 

AM,Q2) = e J ^T^^tn,j,^—^2y-l,jl4y-Aiy^ + j^Aij^l,jl,j^A2j 



The number of gamma matrices to be traced can be reduced with the aid of the identities 

1 , 1 -2p-m^2^± + PlAi,2 
^±^'-'^± {Pi? 

{74, /P-j-} = 2(po + ^JL±ijj) and 74/4i 2 — ~/^i,274- I* follows after some algebra that 

^2 



A4,(gi,Q2) = ^C2(0,0,0;c<.)(qi X q2), 



(B.2) 
(B.3) 

(B.4) 



with 

C2(0,0,0;w) = Svr^ 



(2^) 



P /p Po f 4 2 



(p2)3 (p2)3j (p2)2 (p2)2 



} 



dh 



dz 



(27r)3 In 2Tri 



. 4 



P 



1 



(B.5) 

-(wo -uj)}, 



iQ^/ii — 3 [(z — u;)2 — p^]-^ [{z — uj)^ — p^J^ 

where the angular average of p has been made and the contour C goes around the imaginary axis 



clockwisely. The distribution function /(z) is given by ( 3.13 ). While it is tedious to calculate the 
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residues at the double poles z = ±p and the third order poles z = ±p ± w, there is a short cut 
to extract the infrared divergent piece as cj — 0, that is worth elaborating. Mathematically, the 
contour integral ( B.6 ) should remain finite as a; 0. The divergent terms of the Laurent expansions 
in Lo of all residues should cancel each other. On the other hand, the Baym-Mermin continuation 
amounts to retain the discrete Matsubara lu for the residues so that 



fi±p±Lo)^fi±p) f'i±p±uj)^f{±p) ... (B.6) 



Then the rational dependence left over is continuated to the real w-axis without offsetting (B.6). 
Consequently, the first few terms of the Taylor expansion of f{±p±uj) — /(±p), /'(ipiw) — f'{ip) 
and ... in w, which are required to cancel the small uj divergence are missing, resulting in the infrared 
divergence after the Baym-Mermin continuation. Such terms are easily identified and we find that 

C.(0,0,0,a;) ^—J ^ = - Jjppf {p) = - (B.7) 



which gives rise to ( ^.24 ) 



C. Some technical details behind the one-loop calculation 

In this appendix, we shall expose some technical details behind the one-loop analysis reported in 



the section 3. Substituting eqs.(|3.5|) and (|3.6|) in, the integrand eq.(3.8) may be written as 



^ - (M5 -M5), 



[(Po + - iP' + P5)^ - ^n'^][{Po + - (p - ^5)^ - m"^ 
where p' = p + q and Pg = + 

Itj = -^tr75[-(p' + /i5)74 + z(po + /i - 77174)7 • p']7j[(Po + m)74 - i{p + fJ-bh ' P + n^]7i(C.2) 
+ {P' ^ P;ii^ j) 

= ieijkiPk - P'k)[{p' + M5)(p + M5) + {Po + tJ.){po + m) - m'^] 

and 

J,j = ^tr75[-(p' + ^15)74 + z(po + fi- 77174)7 • p']7j[(po + /^)74 - «(p - ^^5)7 • P + "^]7i (C.3) 
+ (P o P';i^j) 

= ■i<^ijk{Pk+p'k)[{p' + Ai5)(p - /^s) - {p'o + m)(po + ^J■) + m"^]- 
The summation over the Matsubara loop energy pq is straightforward via a contour integral and 
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we obtain that: 



I,j(p,pi5|m) = 



Po 



1^ (C4) 

[{po + w + - (p' + ^5)2 - m^] [{po + - (p + ^5)2 _ m2] 



dz 



2m [(z + a; + A/)2 - Sf] [(^ + _ ^2 j ^,3^ _^ ^ 



7:'^t3k{pk -p'k) f{E'+ - 



(y + /is)^ + ip' + M5)(p + Ms) - 



+ f{E+ - 



jp' + ^5)2 + (p' + /i5)(p + ^5) + ojE'_^ 

E'+{E'+ ~E++ uj){E'^ + E++UJ) 
{P + M5)^ + {p' + /^5)b + M5) + ^^E+ 



E+{E+ -E'++ Lj)(E+ +E'_^+uj) 

jp + Ms)^ + jp' + M5)(P + Ms) - a^£^4 
-E'- uj){E+ + E'-Lj) 



where, the contour C is around the imaginary axis clockwisely. The Baym-Mermin procedure is 
employed to continue the imaginary Matsubara energy uj = 2imTT to the upper edge of the real 
axis to obtain the retarded function. Similarly, 



Ji' 



[(jpQ + cj + /i)2 - [p' + /i5)2 - to2] [(po + - (p - Ms)^ - 



(C.5) 



dz 



c 2m [{Z + LO + fi)^ - Ef] [{z + /x)2 - El] eP- + 1 



-e,-ikiPk+p'k) fiE+ -Ai)- 



-jp' + Ms)^ + jp' + M5)(p + Ms) + 



- /(-i?; - /i) 
+ /(i?_-/i)- 



-(p' + Ms)' + (y + Ms)(p + Ms) - 

-e_+ uj){e'+ + e_+uj) 

(P - Ms)' + (p' + Ms)(p - Ms) - (^E^ 



I{-E+-y) 



E-{E- -E'^+ ^){E- +E'^+uj) 
-(P - Ms)' + (p' + Ms)(p - Ms) + 



E-{E- - E'^- uj){E- +E'^-ij) 



We have then 



n^(Q) 



{Iy(p,^5|m) + Jij{p,fi5\m) ~Itj{p, -Ms|™) - JijiP, -Ms|w) (C.6) 



^ a [Xy(p, Aisles) + Jy(p,Ms|M.) -:J«j(p, -Ms|M.) - Jij(p, -Ms|M,)]}. 



In the static limit, a; = 0, the quantity inside the bracket on RHS of ( |C.5D reduces to R/{p — p') 
and that inside the bracket on RHS of (|C.6D to i?/(p + p') with R the same quantity. We have then 



P-P' , P + P' 



p — p' p + p' p'^ — p^ 



(C.7) 



Upon a shift of the integration momentum, p — > p — q in the terms with E'^ inside the Fermi 
distribution functions, the angular integral becomes elementary. Bearing in mind that Fi{q,Q) is 
real, only the principal part of the angular integral is needed and the result is (3.10) with J-{q\m) 
given by ( 3.11 ). 
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For massless quarks, we may either take the hmit of m — of ( |C.q ) or compute Iij{p, fJ-slO) + 
Jij{p, MsjO) with the massless propagator The contribution from the PV regulators remains 

intact and the result reads 



n^-(Q) 



+ 



PQk + ^Pk 



i2n) 



1 



1 



P 

PV term 



(p - w)2 - (p - q)2 (p + a;)2 - (p + q)2 
1 



/(P-M + Ms) 



{p + w)2 - (p - q)2 {p - cj)2 - (p + q)2 



fi-p + M - Ms) - (M5 ^ -Ms)} 

(C.8) 



with hnuj = where we have shifted the integration momentum according to the prescription 
of the last paragraph. By symmetry, the angular integral may be performed with pk replaced by 
^^Qk and we end up with the form (3.9) with the function Fi{q,uj) given by (3.19). 
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